Evidence for line nodes in the energy gap of the overdoped Ba(Fei_2;Coa;)2As2 from 

low-temperature specific heat measurements 
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Low-temperature specific lieat (SH) is measured on Ba(Fei_a;Coa;)2 As2 single crystals in a wide 
doping region under different magnetic fields. For the overdoped sample, we find the clear evidence 
for the presence of the term in the SH data, which is absent both for the underdoped and optimal 
doped samples, suggesting the presence of line nodes in the energy gap of the overdoped samples. 
Moreover, the field induced electronic specific heat coefficient A'y{H) increases more quickly with the 
field for the overdoped sample than the underdoped and optimal doped ones, giving another support 
to our arguments. Our results suggest that the superconducting gap(s) in the present system may 
have different structures strongly depending on the doping regions. 

PACS numbers: 74.20. Rp, 74.70.Xa, 74.62.Dh, 65.40.Ba 
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I. INTRODUCTION 



Plenty of efforts have been made on the study of the 
iron-pnictide superconductors, since the discovery of su- 
perconductivity with Tc — 26 K in LaFeAsOi-x^xi^ 
An important issue concerning this new family of high- 
Tc superconductors is about the gap structure, which 
should provide clues to the understanding of the micro- 
scopic pairing mechanismj^ Up to date, consensuses have 
been reached on several systems, e.g. LaFePO, KFe2As2, 
BaFe2(Asi_a;Pa:)2, and so on, that nodes exist on the gap 
structurei^^— However, experimental results gave rather 
contradicting conclusions on this issue in other systems 
of the iron-pnictide superconductors One problem 
here may come from the different qualities of the sam- 
ples studied by different groups. Recently, annealing of 
the single crystals in the 122 phase was reported to be an 
effective route to improve the quality of the samples<^"— 
Especially, it was found that annealing can considerably 
decrease the residual SH coefficient 70 in the supercon- 
ducting state and suppress the Schottky anomaly in low 
temperature, which may suggest that fewer impurities 
exist in the annealed samplesi^ This supplies a better 
platform both for investigating the intrinsic properties of 
the samples and for analyzing the SH data. Moreover, 
some groups also point out that the gap structure shows 
a strong evolution with the doping concentrationJ^iiSii^ 

Specific heat is one of the powerful tools to measure the 
quasiparticle density of states (DOS) at the Fermi level, 
so as to detect the information about the gap structure. 
The variation of electronic SH versus temperature and 
magnetic field can be rather different in the supercon- 
ducting materials with different gap structures. From 
textbook knowledge we know that the low-temperature 
electronic SH Cei for a superconductor with an isotropic 
gap should have an exponential temperature dependence, 
namely Cei (x g-^o/ksT^ where Ao is the magnitude 
of the energy gap. However, for a superconductor with 
nodal gap(s), a power law dependence of temperature for 
Cei has been predicted theoretically: Cei oc for the 



gap with line nodes and Cei oc for point nodes i^^'^^ In 
fact, the quadratic term has been observed in cuprate su- 
perconductors, suggesting that line nodes exist in the gap 
functioni2i"i2S In the mixed state, the magnetic vortices in 
superconductors will induce depairing effect within and 
outside the vortex cores, leading to the localized and de- 
localized quasiparticle DOS, respectively. In general, the 
field-induced electronic SH coefficient will show a rapid 
increase with the increase of magnetic field in the low 
field region when the system has the energy gap(s) with 
a rather small minimum value 1^^— 

For the iron-pnictide superconductors, it is rather dif- 
ficult to get the pure electron contribution from the total 
SH, since the upper critical field is very high and the 
normal state can't be obtained by suppressing the super- 
conductivity using a magnetic field. Nevertheless, many 
methods have been used by some groups to subtract 
the phonon contribution from the SH data, which how- 
ever, always inevitably brings about uncertainties J^^i^ 
In this paper, we report the clear evidence of the pres- 
ence of the term in electronic SH of the overdoped 
Ba(Fei_2,Coa;)2As2 single crystals from the raw data, 
which is consistent with the prediction of superconduc- 
tors with line nodes. Moreover, the field induced elec- 
tronic SH coefficient A'j(H) increases very quickly with 
magnetic field for the overdoped sample, suggesting a 
rather small minimum value of the energy gap. This fur- 
ther confirms the conclusion that we stated above. These 
behaviors give a sharp comparison with that observed 
in the underdoped and optimal doped samples, where a 
rather small anisotropy is implied in the gap structure. 



II. EXPERIMENTAL DETAILS AND SAMPLE 
CHARACTERIZATION 



The Ba(Fei_a;Co2;)2As2 single crystals were grown by 
the self-flux method.'^® The as-grown samples were an- 
nealed under high vacuum at 800 °C for 20 days. The 
samples for the SH measurement have typical dimensions 
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FIG. 1: (color online) Temperature dependence of dc magne- 
tization for three Ba(Fei-a;Coa;)2As2 samples with x = 0.065, 
0.08, and 0.13, which are determined to be in the underdoped, 
optimal doped, and overdoped regions, respectively. The data 
are collected with field _ff = 10 Oe using the zero field cooling 
(ZFC) process. The curves are normalized by the magnetiza- 
tion data obtained at 5 K. 
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of 2.5 X 1.5 X 0.2 mm^. The dc magnetization measure- 
ments were done with a superconducting quantum inter- 
ference device (Quantum Design, MPMS7). The specific 
heat were measured with a Helium-3 system based on 
the Quantum Design instrument physical property mea- 
surement system (PPMS). We employed the thermal re- 
laxation technique to perform the specific heat measure- 
ments. The thermometer has been calibrated under dif- 
ferent magnetic fields beforehand. The external magnetic 
field is applied perpendicular to the c axis of the single 
crystals. 

The superconducting transitions of the single crystals 
are checked by the dc magnetization measurements. In 
Fig. 1, we show the temperature dependence of the dc 
magnetization data for Ba(Fei_a:Coa:)2As2 samples with 
nominal doping contents x = 0.065, 0.08 and 0.13. The 
data were collected under a dc field of 10 Oe, which have 
been normalized by the values obtained at 5 K. The sam- 
ple with X — 0.08 is found to be optimally doped with 
the highest onset transition temperature yo"™* « 25.2 
K, which is about 1 K higher than that of the as-grown 
sample. Accordingly, the samples with x = 0.065 and 
0.13 are in the underdoped and overdoped regions, re- 
spectively. 



FIG. 2: (color online) The raw data of the SH for three 
Ba(Fei_iCoa;)2As2 samples with x = 0.065 (a), x = 0.08 (b), 
and a; = 0.13 (c) under different fields in the low temperature 
region. The data are shown in C/T vs plot. We also show 
three solid lines of the theoretical fitting (see text) in each 
figure. 

III. RESULTS AND DISCUSSION 

The raw data of SH for the three samples with x = 
0.065, 0.08, and 0.13 under different fields are plotted as 
C/T vs in Fig. 2. Here we focus on the behaviors 
of our data in the low temperature region below 4.5 K. 
One can see clear different behaviors among the three 
samples. For the underdoped and optimal doped sam- 
ples, the curves display clear linear tendency in this low 
temperature region, which can actually be described by 

C{T,H)^j{H)T + (3T\ (1) 

where 7(i?) is the electronic SH coefficient under a mag- 
netic field of H, and /3 is the phonon SH coefficient. We 
show three typical fitting curves based on equation (1) in 
Fig. 2(a) and (b). From the fitting process, the residual 
SH coefficient under zero field 70 (= 7(0)) is determined 
to be 1.7 mJ/mol and 1.5 niJ/mol for the two 
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samples, respectively. A slight deviation from the linear 
behavior under high magnetic fields in low temperatures 
may come from the influence of the Schottky anomaly of 
the sample platform. We note that the linear behavior is 
rather similar to the previous report. "^^ For the overdoped 
sample, however, the data show a clear negative curva- 
ture rather than a straight line in all the curves up to 9 
T, in sharp contrast with the result of the underdoped 
and optimal doped samples. We argue that this behav- 
ior suggests the presence of the term in the electronic 
SH, which is expected for the superconductors with line 
nodes in the energy gapi ^^i^^ Consequently, the data in 
Fig. 2(c) can be described by 



C(T, H) ^ -i{H)T + a{H)T^ + jST^ 



(2) 



with a{H) the coefficient of the term under the field 
H. As shown by the solid lines, the fitting result is rather 
good. We must note that the negative-curvature fea- 
ture observed here can't be attributed to the Schottky 
anomaly. In general, the behavior of Schottky anomaly in 
SH should display a very strong evolution with magnetic 
field and moreover a steep peak in the low-temperature 
region under low fields should be observed, both of which 
are missing in our experimental data. This gives a direct 
evidence for the presence of line nodes in the energy gap 
of the overdoped Ba(Fei_2;Coa;)2As2. From the fitting, 
we find that the residual term 70 « 1.0 mJ/mol for 
this sample, which is not larger than that of the optimal 
doped sample. This is rather different from the results on 
the as-grown samples.'^'' The situation is reasonable since 
the annealing process can dramatically suppress the num- 
ber of the impurity scatterers in the sample, as pointed by 
K. Gofryk et ali^ Accordingly, we know why the term 
can only be detected in the annealed samples. As shown 
by the theoretical work, the intraband impurity scatter- 
ing tends to make the gap(s) more isotropic, and may 
result in the removal of the low-energy excitations i^i^ 
The fitting parameters for the three samples, /3 and 
a{H), are shown in Fig. 3(a) and (b). One can see 
that the value of (3 is almost independent of field for 
the three samples, which confirms the reliability of our 
analysis and fitting process. Moreover, a{H) of the over- 
doped sample decreases monotonously with the increase 
of magnetic field up to 9 T. For the superconductors 
with a singlet d-wave gap (e.g. the cuprate supercon- 
ductors), it has been pointed out theoretically that the 
term can only exist in moderate temperature regions 
{\/H/Hc2 « T/Tc « 1) in the mixed state because of 
the so-called Volovik effect.— i^i^ Our observation here 
may suggest that some segments of the line nodes are 
not affected by the Volovik effect because the Fermi sur- 
face where they reside in is perpendicular to the magnetic 
field. For the cuprate superconductors, the value of a{H) 
under zero field was estimated to be a(0) ^ ^n/Tc-, with 
7„ the electronic SH coefficient in the normal state.— 
Applying this relation, a(0) is estimated to be about 1 
mJ/mol K^. The value obtained from our data is about 
two times larger than this estimation. This may suggest 
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FIG. 3: (color online) Field dependence of the fitting param- 
eters /3 and ct{H) for the three samples. It is clear that (3 is 
almost independent of the field, while a{H) of the overdoped 
sample decreases monotonously with the increase of the field. 



that the simple relation is not suitable for the present 
multi-gap system, which needs more investigations from 
the theoretical side. In addition, we note that the val- 
ues of /3 decrease with the increase of the doping content 
X. One possible origin of this behavior is that, for the 
underdoped sample, a part of the term comes form 
the contribution of the antiferromagnetic state, which co- 
exists with the superconducting state. There is another 
possibility. As mentioned above, the electronic SH in a 
superconductor with point nodes can also contribute the 
term. So the behavior of the data in Fig. 3(a) may be 
interpreted in terms of the presence of point nodes in the 
energy gap of the underdoped and optimal doped sam- 
ples. However, this argument is rather incompatible with 
the behavior observed in Fig. 4, where the electronic SH 
coefficient increases very slowly with the magnetic field 
both for the underdoped and optimal doped samples (see 
below). One natural argument is that the phonon SH is 
slightly different among the samples with different doping 
contents. 

The field dependence of the electronic SH coefficient 
Aj{H) = j{H) — 7(0) for the three samples is shown in 
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overdoped sample. 

Our results suggest that the gap structure of the 
Ba(Fei_a;Coj;)2As2 system shows a clear evolution with 
the doping content x. Actually, doping dependence of the 
gap structure has been detected by other experiments, in- 
cluding penetration depth, heat transport, specific heat, 
point-contact Andreev reflection, and so onJ^ii^i^id^ 
These results share the same tendency with the present 
work that nodes should exist in the overdoped samples, 
although there is discrepancy concerning the underdoped 
and optimal doped regions. The detailed mechanism of 
this behavior is still vague at this time. Even so, we can 
unambiguously obtain important implications that the 
emergence of line nodes in the overdoped sample is re- 
lated to the evolution of the topology of the Fermi surface 
and the pairing interaction with the doping of electrons 
into the system, rather than the intrinsic symmetry of 
the energy gap. 



FIG. 4: (color online) The magnetic field dependence of 
the low temperature electronic SH coefficient plotted as 
Ay{H)/y„ versus H/Hc2- The values of y„ and Hc2 for 
the three samples are estimated based on the result of other 
groups.— i^i*^ The green dashed line is a guide for eyes, which 
shows the theoretical curve for the samples with an isotropic 
superconducting gap. 

Fig. 4. The data are plotted as Aj{H)/jn versus H/Hc2, 
where Hc2 is the upper critical field. The values of 7„ 
and Hc2 are estimated according to the results of other 
reporta^^i^'^ and shown in this figure. We must point 
out here that the uncertainty of the estimated values of 
7„ and Hc2 will not affect the conclusions deduced in 
the following. In this figure we also display a theoretical 
curve for a conventional superconductor with an isotropic 
superconducting gap by the green dashed line. One can 
see that different behaviors of the three samples in differ- 
ent doping regions are distinct. The data of the optimal 
doped sample seems rather close to the green line under 
low magnetic fields, and even evolves to be below it when 
the field increases. The underdoped sample behaves sim- 
ilarly to the optimal doped sample. On the other hand, 
the data of the overdoped sample increase more quickly 
and are clearly above the green line. The behavior of 
the overdoped sample is quite similar to those reported 
by others, although they have the magnetic field along 
the c-axis of the crystalsf^l We note here that it is not 
easy to describe the data in Fig. 4 using a simple for- 
mula, because the present system was found to display 
multi-band and even multi-gap features. Nevertheless, 
from the present data we can draw the conclusion quali- 
tatively that the energy gap(s) of the overdoped sample 
has rather large anisotropy, compared with that of the 
underdoped and optimal doped samples. Obviously, this 
argument is consistent with our above-mentioned con- 
clusion that line nodes exist in the energy gap of the 



IV. CONCLUDING REMARKS 

In summary, we studied the low-temperature specific 
heat on the Ba(Fei_a;Coa;)2As2 single crystals in a wide 
doping region. Before measurements, the single crystals 
were carefully annealed to improve the sample quality. 
We found a clear evidence for the presence of term 
from the raw SH data on the overdoped sample, which is 
absent both for the underdoped and optimal doped sam- 
ples, suggesting that line nodes must exist in the energy 
gap of the overdoped sample. On the other hand, the 
underdoped and optimal doped samples display rather 
small anisotropy in the gap structure. Moreover, the field 
induced term A'j(H) increased more quickly with mag- 
netic field for the overdoped sample than the underdoped 
and optimal ones, being consistent with the above con- 
clusions. Our results suggest that the superconducting 
gap(s) in the present system may have different struc- 
tures depending on the different doping regions. Future 
investigations on this issue on other systems of the iron- 
pnictide superconductors are needed in the future. 
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